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BESONAiTCE 7IBRATI0NS IB IB" TAKE AND EXHAUST PIPES 



OP IN- LI BE ENGINES* 



Part 3 - The Inlet Process of a Pour-Stroke-Cycle Engine 

By 0. Lutz 



In tho first two parts of this paper (reference l) 
resonance phenomena in pipe lines wore considered* the 
assumption Doing mado that the excitation was producod by 
slmplo volocity fluctuations (which, in tho inlot procoss, 
wore assumod to "be proportional approximately to tho pis- 
ton motion). 

Tho abovo assumption can reprosont only a rough ap- 
proximation which to aome extent enables the resonance 
frequencies to be correctly computed (as in Parts 1 and 
2) but which fails to solve the important problem for the 
inlet process, namely, the effect of the resonance vibra- 
tions on tho volumetric efficiency. If tho lattor is af- 
foctod by tho vibrations, the volocity fluctuations at 
tho point of excitation aro no longor proportional to tho 
piston motion but aro a function of tho vibrations. It 
is evident that any mode of computation which does not 
tako this into account can not bo used to oxplain this 
phonononon, although previously this had always boon at- 
tempted. 

Taking tho case of rosonance vibrations which are 
excited by a periodically controlled opening, tho author 
had developed a method (roferonco 2) in which tho boundary 
conditions woro made to correspond moro noarly to actual 
conditions than has been dono hcretoforo and by which the 
change in tho flow due to tho fluctuations could be com- 
puted. In the present paper this method will be extended 
to the inlet process of a four-stroke-cycle engine. In 
section 1 the previously developed method will be briefly 
described. In section 2 application will be made to the 



*Uber Besonansschwingungen in don Ansaug- und Auspuff- 
leitungen von Boihenmotoron. 3. Teilbericht: Bio Anregungs 
bedingungon beim Ansaugvorgang von Vior taktmotoren. Luft- 
f ahrtf orschung, vol. 17, no. 4, April 20, 1940, pp. 123-28. 
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inlot process and the required "boundary conditions ob- 
tained. In section 3 the conditions obtained will he 
applied to the simplest case of a single-cylinder engine 
and, finally, in section 4 the basic principles will be 
given for application to multicylinder arrangements. 



1. PEHTCIFLE OP THE METHOD 

The case considered is that sketched in figure 1. 
At the end of the pipe x = I, the flow passes through 
a periodically controlled opening into a chamber whero no 
back pressure is assumod to be built up. Tho narrowost 
flow cross soction is donotod by f, the pipe cross soc- 
tion by f r . With the notation of figure 1 we have 

|(v 8 - r r ») = Pr - p 

(whore the local acceloration term is neglected). Using 
tho continuity condition: 

f v = f r v r 

wo havo 

v r = ■■ f /-(p r - p) (1) 

Defining nean values v , f Q , p r0 , and p 0 by 

and noglocti:!;-; products and powers of the rolativo devia- 
tions from tho mean valuos etc., wo have 

v 0 f o 
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and 



Av = — » — (1 + — r— * r-. . .) + v- ( —7— 4 r. . .) 

f o ^ 2(p ro - p 0 ) / 0 \2ivro - Po) * 

(lor simplification J f r 3 - f 3 has "been set equal to 



f r a - f c a Bih.ce the Jet cross section f will always 



J 

be small compared to the pipe cross Beotlon f r .) 

In the papor Cited previously the idea was developed that 
for tho computation of tho magnitude and phaso of tho vi- 
bration the linear terms are used while for tho computa- 
tion of tho changes in tho moan va.luos tho torms of highor 
order woro used, terms of seoond order only being taken 
into account at first. Thus, setting_the mean velocity 
equal to the sum of the mean value v Q from equation (s) 
for the system without vibrations and- the increment Av 0 , 
which gives the changes by the vibration process: 



v Q .= v Q + Av 0 (4) 



we have 



v o 



2fQ^Pro ~ Po^ 

As described in the first part of this paper (reference l) 
and in the previous work, tho computation is carried out. 
using complex variables, sotting 



Av = E V 



incut 



Af = £°- S j*e lauJt 
Ap r = pa S P ?n e lnUJt 
Ap = pa S.£*o illU,t j 



(7) 



(a denotes the velocity of sound; all vibration vectors 
are converted to the dimensions of a velocity.) 
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Iron equation (5) there follows the "boundary condition 

In = 3£ +£ (Z*n " W (8) 
where € is the "opening coef f ioionf n given "by 

g= 9 a " ?0 - , " ^ (9) 

2(? ro - P 0 ) ^ 0 ( f r a " f o a ) 

Tor the deviation av 0 from the mean value there is ob- 
tained from equation (6) 

j. _ 2TT , 



o = tct 

£ 1 ^a^^M, - p*)o in(ut l dt do) 



Av o v 0 t 0 / V***" /' : *• k -rn in'- j 
In tho previous work it is then provon that 

Av 0 = ZP*(P? a - W (10a) 

* T o 

where the sum is "built up of scalar products. Thus, if 
the difference vector Z*n~ — n is ou * of P^ase "by less 
than a quarter period with the opening vector JjJ, the 
average quantity of charge delivered is increased by the 
vibration process. In the previous work examples for 
this case of excitation were computed and agreement of 
tho toot rosults with the computation proven. 

2. APPLICATION 03" THE METHOD TO THE IHLET PEOCESS 



Tho inlet procoss in tho caso of a f our-s troko-cyelo 
engine can be divided up into two proeossos, tho first of 
which was alroady considorod in soction 1, namely, tho 
flow process through tho controllod opening and tho charg- 
ing process inside tho cylinder. 

Tor the first procoss the boundary condition is al- 
ready known: 

I r = 1* + € (P* - P*) (8a) 

(written down for the fundamental); tho socond procoss 
requires furthor consideration. 
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Tor an adiabatic charging process, we' hare, with the 
-notation, of figure 2X. 

to b k *S . K M. (n) 

where G is the variable weight of the gas t J the in- 
stantaneous cylinder volume. - The rate of change with time 
is 

dt dt Kp dt 

Converting- the cylinder volume into pipo lengths by 



there is obtained with = Yf r v,.: 

dt r r 

v r = 4£-+JL4f (13) 
r dt up o.t 

If no pressure fluctuations were set up in tho cylindor, 
the boundary condition would be 



T r = 



ds 
dt 



(12a) 



As In section l r the mean values and the deviations from 
those moan valuos 



ds 
dt 



.c 0 (l + |^) .-..(1**J) 



are again substitutod in equation (12) 

As 

1 + — 

b 0 d A p 

T o f ' " " ' 



(i + *SN = 0o (i + + JL°. ?J> £Ap 

\ *V °\ c 0( / KP 0 j_ h Ap dt 



and tho vibration magnitudos as in oquation (7) expressed 
in complex form: 
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* • r ■ r ■ " - 1 ' ■ . 

■ J .... 1 

■ ■ 1 



Ac = EC^o 1 ^* 
As =Zsle liwt 



(13) 



The linear tores thon lead to the 'boundary condition 
(written for tho nth harmonic) : 



T„ u = C„ + i 2. p« 



(14) 



The torus of second order give for tho moan values: 

iML "t Ap ^ 
dt - P aS P 0 



v o ■ c o + T7a As — »„a^ a P 7f (15) 



In the same manner as for equations (10) and (10a) , setting 
Av Q = v 0 - c 0 , there is finally obtained 

Av 0 . i n P* (15a) 

(The last term in equation (15) gives no contribution 
since the first derivative of the function itself leads 
by a quarter period, so that the scalar product of the 
two vanishes.) 

Betveen the harmonics 5^ , vrhich give the cylinder 
volume at any instant, and tho harmonics 0^ , which give 
the instantaneous piston velocity, thero is a definite ro- 
lation. In figuro 5 tho periodic relation is drawn for a 
f our-stroko-cycle onginc (effect of finito connecting rod 
noglocted). Tho fundamental frequoncy of tho inlot proc- 
oss (that is, half tho ongine froqxiency) is donotod by (0 . 

We have 

o = 77 , thus c x = 2uj 8 1 (16) 
dt 

Expansion into a harmonic series gives 
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a =' b' + ^ s i sinuj t + ~ a'in'^u) t + a x Bin 3uj t 

0 3tt 4 .5tt 



+ — s, sin 4u> t 
3tt 1 



and 
c 



= ~ c. + c, cbato t + cob 2u) t 

2tt 1 3tt 1 4 



a/2 „ 1 . , 

+ -—- Ci cos 3 cut + — o, cob 4 U)t...« 
5tt 3tt 



(17) 



from which, using equation (16), there is obtained 



dS n n a 

■■ = — G 
dt 4 



(18) 



or 

Equation (15a) then goes over into 



(18a) 



(15b) 



The boundary condition at the cylinder end of the pipe 
is thuo given by the vector equations 



Irn 



25 + e <P$ a - P*) 



and 



— rn - 



= C + i p* 



(8) 



(14) 



while the chango in the volumetric efficiency due to the 
vibration procoss is given by 

.•a 



av = - _ v C P* 
o Qa ^— n — n 



(15b) 



It should be pointed out once more that n. gives the order 
of the harmonic and is not a rotational speod. 



3otv/oon tho "opening coefficient" € in equation (8) 
and the frequency uj for the inlet process of a four- 
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stroke-cycle engine there i-e,_a hyperbolic .relation. For, 
from equation (9), 

a f 0 3 
▼ 0 (fp" '-"*o") 

The moan pipe velocity' 1 v (without taking account of the 
vibration effect) is proportional, however, to the engine 
3pcod and thua to the fundamental frequency u) • Thus, 
from equation (17), 

whonco, using oquation (16), 

_ uusi 
v o = — 

According tc figure 3, the stroko amplitude s and tho 
moan valuo s 0 arc connected by the relation 

s 0 + b x 

E - R = £ V 

( € v is tho compression ratio), so that, finally, 

€ w - 1 U)S_ 

and for tho opening coefficient thoro is obtained 

<=-5§7 tl9) 

a 

with 



= _^(lv + Dfo a (aQ) 

<C„ - l)(f » - f ■) 



3. EXAKFLJE:. S I IT GLE- CYLINDER IHGIKE 
Lot tho pipo connocting to the slnglo-cylindor engine 
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oe of length I (fig. 4). As proven in the first paper, 
for. an open pipe the following relations exist "between 
the vibration magnitudes 



and 



Io 



i +«« "Li v 
- i tan — I r 



COB 



(21) 



tfith "boundary conditions (8) and (14) , written down for 
the fundamental, and taking account of (19), there is ob- 
tained after elimination of the magnitudes PS, P* , and 



Ir* 



Io (cos *± - ^ B m i£ + ^(212. y cos i£L ) 

m Q + i(™J)\* (22) 
B \ a /— 

for s Q = 0, the equation assumes the known form 



Io C09^= 0 (22a) 

Equation (22) describes the vibratijon processes in a 
general form. If the dimensions of the system are known, 
the entire process can "be followed with tho aid of this 
equation. (In a lator report, an example will ho computed.) 

The conclusions in the prosent paper .will be obtained 
witli the aid of a simplified example. 

It will bo assumed that 

s 0 = and " p = 1 (23) 

and, further, that the controlled opening follows the vari- 
ation in the piston speed in figure 3. Since all the vi- 
bration magnitudes are roduced to the dimensions of a ve- 
locity, we have in this caee simply 

' 1* ■ C. (23a) 
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Finally, let the nondlmens ional expression — = cp. The 

"boundary conditions (8) and (14) then become (for the 
fundamental) 

and 

I- = C | P* (14a) 

and ve have in place of (22) 

X 0 (cos cp - | sin cp + i ^-cos - J* (i + 1 ) <22b) 

from which, with the aid of (21), there 1b obtained 

1 - ftancp + SL. *E? tancp 

I r " rC-i- 2 f 0 (24) 

fl-2.tan (l -9£tanq)) 3 + f- 

\ 2 > 16 ^2 / 16 

From (14a) there ia further obtained 
a 

- SL_ tan cp tan cp - T. tan B tp 

Z* = 7 ra s- 0 - i ^ ^-r 5- C (25) 

(l-fto«p) + SL (! - | tOB „)" + 2- 

From the above equation (25) and equation- (15b) tho effect 
of tho fundamental on the change in tho volumetric effi- ■ 
cioncy can bo obtained. 

Wo havo 

CP* 

Av ° = - ST 

which, with the aid of equation (25) bocomos 

32a" tan * „ a 
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where C x denotes the absolute value of tho amplitude of 
the fundamental vibration* 



Trom equation (17) 



and from (17a) 



Oi = -Z— ci 
3tt 



4/2 - 
1 C 3 V ° 



and using tho relations that led to formula (19) 

"~ v ^ 8 (26) 



v 0 9tt £ v + 1 



"he rosults can host ho represented graphically. In fig- 
uro 5 the relation botweon the piston velocity fluctuation 
and tho volocity fluctuation V r according to (24) is 
givon by a polar diagram. Tho deviation from tho usual 
assumption that 7 r is oqual to C is, for tho oxaraplo 
chosen, often cons idorabl o . Tho phaso shifts can amount 
to moro than a quarter period and considerable difforoncos 
arise also in tho omplitudos which difforoncos may bo par- 
ticularly largo in tho froquoncy rango boforo tho first 
resonance (0 < rp < 90°). 

Jiguro 6 shows tho variation in tho velocity fluc- 
tuation 7 0 at tho pipe ontranco (according to (21)). 
It is interesting that no narkod rosonancos occur although 
no damping has boon assumed. Who re as , with tho usual as- 
sumption 7 r = S, the velocity fluctuation V Q 1 is in 
phase with C_ Xshown dotted) and becomes infinitely great 
at tho rosonanco positions (90°, 270°) tho phase position 
of T 0 changes ■ continuously but no infinitely largo ampli- 
tudes any longor occur. The oval-shapod loci describod 
by the tip of tho veotor J Q toko a surprising course. 
The first oval, which is for 0 < cp < 180°, is largo; tho 
socond (180° < «p < 360°) is snail. Iron then on, how- 
over, tho ovals bocomc increasingly largo. Examining tho 
relations more closoly, wo find that the valve motion is 
mainly responsible for this Increase in tho vibration val- 
ues. It may thoreforo bo concludod that oach change in 
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the valve notion oust have an inportant effect on the 
process. lor largo f requencies , that is, for large val- 
ues of cp or U)l/a, the linear and tho quadratic terns 
in oquation (22) are predominant, so that oquation (22) 
nay "bo replaced "by 



or 



\ a u js 0 a / 

a 

Without giving a detailed proof (reference 2 , p. 37), an 
equation of the forn (22c) describes the process at the 
entrance of a throttled pipe line (fif. 7) which at the 
other end is excited by a fluctuation in the size of the 

g 

opening j?*. The factor — — corrosponds to the"throt- 
~ ujs 0 

~a 

tling coofficiont" 8 X , which, in the refcrcnco cited, 
was derived as 

o £r!lo_ V 
1 ~ f^ a av 0 

At the higher frequencies therefore the system may be re- 
placed "by that shown in figure 7, where the throttling, it 
may he remarked, is to bo chosen smaller with increasing 
frequencies. 

In figure 6, finally, the variation of tho volumetric 

av Q _ U) I 

efficiency -~ — is. shown as a function of <p » -jj-, if a 

To" 

compression ratio of C T - 6 is assumed* It is found 
that the vibration may t,-ivo riso to improvement as well 
as to impairment of .the volumentric efficiency. The am- 
plitudes _V 0 of tho velocity fluctuations at the pipe 
ontranco (from fig. 6) aro also drawn. Tho first roso- 
nanco position occurs ahead of the valuo (90°) valid for 
tho assumption T r = G_, the difforenco amounting to 33 
pcrcont. The sub~noquont resonance positions, however, 
diffor only slightly (270°, 450°,' 6300...). Tho figuro 
also shows cloarly the increase in tho. maximum valuo at 
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the highor resonance positions (lesser throttling). The 
gains in volumetric efficiency always occur ahead of the 
vibration resonance positions, "become equal to zero at the 
resonance positions, and bocome negative after the rosonanco 
positions. In the oxamplo oheson, the host value of the 
gain is about 30 percent and at the highor rosonanco posi- 
tions no ffloro than 6 or 7 percent. If the fact is taken 
into account that the vibration process is assumed to bo 
without lossos (f rictional and damping) , no important of- 
f octs should bo expected, at least at the high resonances. 

A later report will give the results of an investiga- 
tion made with the aid of examples to determine how strong 
is the effect of certain factors, particularly a change in 
the timing of the opening cross section. 

B"or the present the following conclusions may be 
drawn: 

1. Tho inlet process in the caso of a four-stroke- 

cycle engino loads to vibrations which may 
considerably diffor from thoso that are ox- 
citod by a puro volocity fluctuation. With 
opon pipes no infinitely largo vibrations occur. 

2. Tho froquoncios of tho largo vibrations (resonance 

froquoncios) are lowor at tho first rosonanco 
position than in the case of volocity excitation 
and tho intake pipe thus appears as olongatod. 
The higher rosonances, howovor, are in good 
agroomont. 

3. At highor froquoncios tho systom operates liko a 

pipe throttled at tho ontranco, which at its 
other ond is oxcitod with a velocity fluctuation- 
equal to tho opening f luctuat ion, and it is neces- 
sary that tho throttling be reduced with increas- 
ing froquoncy. 

4. As a rosult of the vibration the volumetric offi- 

cioncy can bo considerably affected. 

5. A gain in tho volumetric efficiency always tokos 

placo ahoad of tho rosonanco positions (in the 
sonso of increasing values of (ul/a). It is 
largest boforo tho first resonance position (30 
porcont for tho example chosen). Beyond the 
resonance positions tho volumetric efficioncy is 
impaired. 
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4. MULT ICTIIH33EE EUG-IITES 



The treatment of multicylinder engines presents no 
fundamental difficulties. The phenomena soon "become so 
complicated, however, that the general treatment promises 
few results of fundamental importance. The individual 
treatment of a given system is not difficult. The rela- 
tions for a two-cylinder system will ho briefly devolopod. 

The intake piping of a two-cylindor engine is 
skctchod in figure 9; It is known from the previous pa- 
per that the process can bo simply followed if tho pipes 
are considered as made up of simple "open" pipe parts, as 
shown in figure 10. Pipe 1 (length I ) has the required 
velocity fluctuation 7 0 at its entrance; the pressure 
fluctuation P Q * can be set equal to zero (open pipe). 
At the end of this pipe the fluctuations are given by the 
following known relations 



P* « 

£r3 



= - i 



sin Bli T 0 
a 



-II 



cos ^ T 0 
a — 0 



Pipe 2 (loagth 1 3 ) is also "open" at tho entrance since 
at this point the pressure fluctuation to be superlmposod 
on that of pipe 1 is PJ X = 0; the exciting velocity fluc- 
tuation (to bo determined) is T_ . Corresponding to thoBO 
initial conditions , the fluctuations at tho other end are 
given by 

p»i = + i ain v . 

— r 3 a * 

V n = - cos Yj 

Tho fluctuations of the ontiro system at tho other end 
aro theroforo 



(27a) 

III = cos ~- lo ~ cos V T (27b) 
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It 1b now necossary to find the relation betweon the flue* 
tuation values of the velocity, tho proBsmo, and the in- 
tako system (piston velocity and control opening). . 

Boundary conditions (8) and (14) using the opening 
coefficient e given by equation (19) , may bo written 
(for tho fundamental): 

1 T = 1* + — ^— (Zr* - 2*) (8b) 



and 



IDS 

a 



IDs _ 

V_ = 0 + i 2. P* (14b) 

A ft. 



Eliminating P * , tho rolation botweon V r and „P r *» 

C + i ° 1* i 

Z r &5 + 2 P r * (28) 

1 + i ^ 1 + i 

Pa a Pa 3 

and sotting 

C + i =— P* 

Pa 3 " 

& 2 = (39) 

1 + i 



and ■ 



Pa s ■ 



* = V"7 ■ (29a) 

1 + i^ S ° 



Pa a 

tho rolation may bo put in tho gonoral form 



It = £ £ + * P * (28a) 



If oqual cylinders aro assumed, with tho arrangement 
of figure 9, the values a and b for both pointB of 
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oxcitation aro equal, and the following relations aro ob- 
tained 

Vj = a Cj + b P^ (27c) 

and 

V ZI = a C IZ + b P* 3 (27d) 

Tho pressure fluctuation which occurs in the open pipe 
portion 1 at the distance l x from the pipe entrance is 
given "by 

Ux " - 1 8ln ^ *o < 27e > 

From the five relations ((27a) to (27e)) thus set up, the 
five unknowns J_ 0> V It V IIt P£ lt and £ r * may "be deter- 
mined. 

In the same manner every system aan "be investigated. 
General conclusions are difficult to draw because the co- 
efficients a and b are directed magnitudes (vectors) 
which are not constant but vary with the frequency. 



5. SUMMARY 



Using a previously developed method, the "boundary 
conditions for tho resonance vibrations during the inlet 
process of f our-stroko-cycle engines aro set up. Tho 
results deviate considerably from those obtained under 
the assumption that the velocity fluctuation is propor- 
tional to the cylinder piston motion. Tho deviation is 
less at tho position of resonanco froquencics. By tho 
method developod, tho offoct of the resonance vibrations 
on tho volumetric efficiency can be demonstrated. 



Translation by S.. Roiss, 
national Advisory Comnitteo 
for Aoronautics. 
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Figure 2.- Sketch 
for 

Betting up the 
boundary condition 
for the inlet 
procosB. 
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Figure l.~ Sketch for 
setting up 
the boundary condition 
(8). 
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c=Cj_ cob 2wt b=6q+ Bj_ sin 2wt 
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Figure 3.- Piston velocity and cylinder volume during the inlet process. 




Figure 4.- Intake pipe 
to a single- 
cylinder engine. 
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Figure 6.- Velocity fluctuation V 0 at 

the pipe entrance (V 0 * holds 
Figure 7.- Equivalent system. for the usual assumption V r = C) . 
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Figure 10*- Uethod of considering 
the pipe in Fig. 9 as 
divided into two open pipe parts. 
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